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Lateral nozzle forces are known to cause severe structural damage to any new rocket 
engine in development. Currently there is no fully coupled computational tool to analyze this 
fluid/structure interaction process. The objective of this study was to develop a fully coupled 
aeroelastic modeling capability to describe the fluid/structure interaction process during the 
transient nozzle operations. The aeroelastic model composes of three components: the 
computational fluid dynamics component based on an unstructured-grid, pressure-based 
computational fluid dynamics formulation, the computational structural dynamics 
component developed in the framework of modal analysis, and the fluid-structural interface 
component. The developed aeroelastic model was applied to the transient nozzle startup 
process of the Space Shuttle Main Engine at sea level. The computed nozzle side loads and 
the axial nozzle wall pressure profiles from the aeroelastic nozzle are compared with those of 
the published rigid nozzle results, and the impact of the fluid/structure interaction on nozzle 
side loads is interrogated and presented. 
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Nomenclature 


Ci,C 2 ,Cs,C f i= turbulence modeling constants, 1.15, 1.9, 0.25, and 0.09. 

C p = heat capacity 

D = diffusivity 

F yZi F y> F z = integrated force, and component forces in the lateral direction 

/ = frequency 

FI = total enthalpy 

K = thermal conductivity 

k = turbulent kinetic energy 

L/S = ratio of long-axis to short-axis 

Q = heat flux 

T = temperature 

t = time, s 

u = mean velocities 

V 2 = Z u 2 

x = Cartesian coordinates or nondimensional distance 

a = species mass fraction 

e = turbulent kinetic energy dissipation rate 

6 = energy dissipation contribution 

= viscosity 

= turbulent eddy viscosity (=pC M k 2 /s) 

FI = turbulent kinetic energy production 

p = density 

a = turbulence modeling constants, 0.9, 0.9, 0.89, and 1.15 for Eqs. (2), (4), (5), (6). 

r = shear stress 

co = chemical species production rate 

Subscripts 
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r 


= radiation 


t = turbulent flow 


I. Introduction 

Nozzle lateral forces during engine startup and shutdown transients, if not properly managed, are known to cause 
severe structural damage to the engine hardware for almost all liquid rocket engines during their initial development 
[1-4]. Transient nozzle side load is therefore considered a high risk item and a critical design issue. For that reason, 
many research efforts [5-19] have been devoted to understanding the side load physics and their impact on the 
magnitude of side loads. For regeneratively-cooled engines such as the Space Shuttle Main Engine (SSME), the 
peak side load generating physics have been identified as the X shock oscillation across the nozzle lip [7]. For film- 
cooled engines such as the Japanese LE-7A engine and the U.S. J-2X engine, the major side load generating physics 
have been associated with the jump of the separation line [3, 8]. Other side load physics such as the Free-Shock 
Separation (FSS)-to-Restricted-Shock Separation (RSS) transition have been mentioned as the critical physics for 
the European Vulcain engine [9], 

In the aforementioned research efforts, computational fluid dynamics (CFD) have been demonstrated as a 
powerful analysis and design tool in computing and understanding the underlying transient side load physics inside a 
rigid nozzle. During an actual hot-firing of an engine, the nozzle wall or the structure of the nozzle, flexes or 
deforms in response to the lateral aerodynamic forces. The deformation of the nozzle wall simultaneously modifies 
the aerodynamic flowfield and the lateral forces, which in turn affects the nozzle wall deformation. This aeroelastic 
movement of the nozzle wall, which was not considered in the rigid nozzle modeling, is one of the important side 
load physics and needs to be considered. 

The aeroelastic modeling of the transient nozzle side loads, requires the coupling of the CFD and computational 
structure dynamics (CSD) in order to describe the fluid -structure interactions (FSI). To the best of the authors’ 
knowledge, there are only two recent exploratory papers [20, 21] that attempted to address this issue. Ref. 20 
demonstrated the displacement of a two-dimensional, smooth nozzle wall of the J-2S engine at a fixed pressure ratio. 
Ref. 21 performed a quasi-static FSI simulation of a three-dimensional, smooth nozzle wall of the SSME. The 
computed elapse time, however, was only 0.02 s (out of a 5 s nominal startup transient) and the timing of the 
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simulation was away from the two major side load physics, FSS-to-RSS transition and X shock oscillation across the 
nozzle lip. 

In this effort and our first attempt on this subject, an aeroelastic model is being developed that composes of three 
components: the CFD component based on an unstructured-grid, pressure-based formulation, the CSD component 
developed in the framework of modal analysis, and the fluid-structural interface component. A three-dimensional 
simulation of the transient startup sequence of the SSME is being performed. The computed results with this 
aeroelastic model will be compared with those of the published SSME rigid nozzle effort [7], 


II. Computational Methodology 

A. Computational Fluid Dynamics 

The computational fluid dynamics (CFD) methodology is based on a multi-dimensional, finite-volume, viscous, 
chemically reacting, unstructured grid, and pressure-based formulation. Time-varying transport equations of 
continuity, species continuity , momentum, total enthalpy, turbulent kinetic energy, and turbulent kinetic energy 
dissipation were solved using a time-marching sub-iteration scheme and are written as: 
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A predictor and corrector solution algorithm was employed to provide coupling of the governing equations. A 
second-order central-difference scheme was employed to discretize the diffusion fluxes and source terms. For the 
convective terms, a second-order upwind total variation diminishing difference scheme was used. To enhance the 
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temporal accuracy, a second-order backward difference scheme was employed to discretize the temporal terms. 
Point-implicit method was used to solve the chemical species source terms. Sub-iterations within a time step were 
used for driving the system of second-order time-accurate equations to convergence. Details of the numerical 
algorithm can be found in Refs [22-25]. 

An extended k-s turbulence model [26] was used to describe the turbulence. A modified wall function approach 
was employed to provide wall boundary layer solutions that are less sensitive to the near-wall grid spacing. 
Consequently, the model has combined the advantages of both the integrated-to-the-wall approach and the 
conventional law-of-the-wall approach by incorporating a complete velocity profile and a universal temperature 
profile [27]. A 7-species, 9-reaction detailed mechanism [27] was used to describe the finite-rate, hydrogen/oxygen 
afterburning combustion kinetics. The seven species are H 2 , 0 2 , H 2 0, O, H, OH, and N 2 . The thermodynamic 
properties of the individual species are functions of temperature. The multiphysics pertinent to this study have been 
anchored in earlier efforts, e.g., SSME axial force and wall heat transfer [22], SSME startup side load and dominant 
shock breathing frequency [7], J-2X startup and shutdown side loads for a nozzlette configuration [8], nozzle film 
cooling applications [28], and conjugate heat transfer [29], 

B. Computational Structural Dynamics 

The structural dynamics response due to fluid flow actions has been analyzed using direct finite- 
element analysis. The aeroelastic equations of motion of the solid bodies are given by 


[m]{f}+[c]{f}+[x]{f}={f} (7) 

where {F} is the displacement vector, [M] is the mass matrix, [C] is the damping matrix, [X] is the 
stiffness matrix, and {F} is the force vector due to the aerodynamic loads and shear stresses. 

The motion Equation (7) of the structure can be solved using modal approach. On the basis of modal 
decomposition of the structure motion with the eigenvector of the vibration problem, the displacement, 
velocity and acceleration can be transformed to the generalized displacement, velocity and acceleration 
using a transformation matrix, which can expressed as the following: 
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jy}= [®]{z} ; {rj= [®]{z}. {r}= [®]{zj 


( 8 ) 


Here [®] is the mode shape matrix containing the eigenvectors, orthonormalized with the mass matrix, 
{z}, jzj, and jz} are the generalized displacement, velocity and acceleration vectors, respectively. The 
eigenvectors are orthogonal to both mass and stiffness matrixes and if Rayleigh damping is assumed, it is 
also orthogonal to the damping matrix. Pre-multiplying Equation (8) by [®] , we get 


{z}+ [®f [c][®]{z}+ [®f [k]Mz} = [of {F} 


(9) 


where 




Equation (3) can be written as n individual equations, one for each mode, as follows 


fzi+2g®,z, +a>?z, =r t 


i = l,2,...,n 


( 10 ) 


Here C 0 j is the natural frequency for the ith mode c, is the corresponding damping parameter for that 
mode. The solution to Equation (10) can be obtained for each mode using direct integration algorithm. 

C. Transient Startup Sequences 

Transient system-level simulation is a vital part of the computational methodology, because it provides the time- 
histories of the inflow properties entering the nozzle. Simply put, the ramp rates, or histories, of the inlet pressure. 
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fluid temperature and species concentrations play an important role in determining the type of side load physics, 
magnitude and duration of the side loads during the transient operations. In other words, the time-varying inlet flow 
properties determine the residence times of the side load generating flow physics inside the nozzle. 

The system-level simulation is based on a lumped, control-volume approach to model the rocket engine as a 
network of components and sub-components. This method of transient system-level modeling has been shown to be 
effective in simulating the low-frequency, transient physics associated with the operation of previous and existing 
rocket engines (SSME, RL-10, IPD, etc.) and therefore, is an important tool in the design and planning of 
sequencing the transient events of rocket engine operation. Figure 1 shows some of the inlet flow properties 
obtained from the system model: the time-varying inlet pressure, temperature, and equivalence ratio profiles. These 
time-varying inlet properties were used at the injector faceplate of the thrust chamber for the CFD computation. Two 
significant pressure rise events can be identified in the inlet pressure history of Fig. 1. The first one occurs at 1.5 s 
due to oxygen prime, while the second one occurs 
at about 2.4 s, caused by the step opening of the 
oxygen valves in the pre -burners. The inlet 
temperature history shows a sharp jump at 1.5 s, 
leveling off after 1.75 s, jumps a little bit again at 
2.4s, and increases linearly until around 3.1 s when it 
reaches the final temperature. The inlet equivalence 
ratio history shows that the thruster environment is 
fuel rich throughout the start-up transient, especially 
in the first 1.5 s, setting up the potential for 
afterburning. That turns out to be the source of the 
combustion wave, because the pressure jump at 1.5 s 
increases the reaction rate of afterburning, which 
leads to the generation of the combustion wave. Afterburning plays an important part in the subsequent asymmetric flow 
physics such as the shock transitions and shock breathings across the nozzle lip. As mentioned in the beginning of this 
section, that the route or history between the starting and end points of any of the curves in Fig. 1 influences the side load 
physics intimately, any simplification on any part of the sequence may run the risk of missing or degradation of 



Fig. 1 Computed thruster chamber inlet properties during 
the start-up transient. 
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important side load physics. This SSME start transient process involves thermal-fluid physics phenomena and safety- 


based operating practices that are typical of a conventional LH2/LOX rocket engine. 


III. Computational Grid Generation 

This section will be completed when the final manuscript is due. 


IV. Boundary and Inlet Conditions 

Since SSME is a first stage engine, fixed freestream boundary conditions were set corresponding to sea level. 
Time-varying inlet flow boundary conditions were used at the main combustion chamber (MCC) inlet. These inlet 
flow properties (obtained from the engine system model) were the time varying total pressure, temperature and 
propellant compositions. For engine startup computations, the thermal wall boundary condition was initially 
adiabatic. When the startup transient simulation reached 1 .5 s, the point where the MCC pressure started to ramp up, 
wall temperature profiles obtained from a separate steady-state calculation were imposed onto the MCC and 
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Fig. 2 Computed side forces for the rigid nozzle case. 


regeneratively cooled nozzle wall. 


combustion wave 


lip RS: 


Mach disk flow: 
FSS 


core jet flow 




FSS-to-RSS 

I 

Mach disk flow: RSS 




S oscillation 


nozzle flowing full 


V. Preliminary Results and Discussion 

The computations were performed on a cluster machine using 10 to 15 processors. For these transient 


computations at sea level, the reaction rates of combustion are high and global time steps ranging from 1 to 10 ps 




are expected to be used throughout the computations. These global time steps used correspond to CFL numbers 
ranging approximately from 0. 1 to unity. 

The Rigid Nozzle Case 

Figure 2 shows the computed time-varying side loads and associated physics, simulating the SSME rigid nozzle 
hot-fired at sea level [7]. The side load physics captured included the combustion wave, FSS-to-RSS transition, and 
lip RSS oscillation (or breathing). 



Fig. 3 Computed first three nodal diameter modes. 


The Aeroelastic Nozzle Case 

The aeroelastic-body computation is currently in progress. Figure 3 shows the deformation of the SSME nozzle 
due to the first-, second-, and third-nodal diameter modes. It is expected that the first-diameter mode, or the 
ovalization of the nozzle dominates the fluid-structure interaction during the transient startup process. At present 
moment, the modeling of the frequencies of the first three nodal diameter modes has been completed and a 
preliminary computation of the transient startup process is in progress. 

VI. Conclusion 

The conclusion of the results will be ready when the final manuscript is due. 
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